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We show that spin squeezing implies entanglement for quantum tripartite-state, where the sub- 
system of the bipartite-state is identical. We study the relation between spin squeezing parameters 
and entanglement through the quantum entropy of a system starts initially in a pure state when 
the cavity is binomial. We show that spin squeezing can be a convenient tool to give some insight 
into the subsystems entanglement dynamics when the bipartite subsystem interacts simultaneously 
with the cavity field subsystem, specially when the interaction occurs off-resonantly without and 
with a nonlinear medium contained in the cavity field subsystem. We illustrate that, in case of large 
off-resonance interaction, spin squeezing clarifies the properties of entanglement almost with full 
success. However, it is not a general rule when the cavity is assumed to be filled with a non-linear 
medium. In this case, we illustrate that the insight into entanglement dynamics becomes more 
5^ ' clearly in case of a weak nonlinear medium than in strong nonlinear medium. In parallel, the role 

of the phase space distribution in quantifying entanglement is also studied. The numerical results 
of Husimi Q-function show that the integer strength of the nonlinear medium produces Schrodinger 
cat states which is necessary for quantum entanglement. 
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I. OVERVIEW 

The squeezing of the quantum fluctuations is one of the most fundamental manifestations of the Heisenberg uncer- 
O" 1 tainty relation, which is among the most important principles of quantum mechanics. A long time ago, great effort 
has been paid to squeezing of the radiation field due to its strong application in the optical communication [l| and 
£«C) . weak signal detections [2|. Accordingly, it was established the relationship between the squeezing of the atoms and 
J> ' that of the radiation field Q and the possibility of squeezed atoms to radiate a squeezed field [jj. Moreover, much 



attention has been devoted to atomic spin squeezing |5l— fl2jj - Spin squeezed states (ci 19.113b 27] are quantum correlated 
states with reduced fluctuations in one of the collective components. Sp in squeezed states offer an interesting possi- 
bility for reducing quantum noise in precision measurements [H, [3 [28^3^ , with potentially possible applications in 
NJ . • ultra- precise spectroscopy, atomic interferometers and high precision atomic clocks [3 . 

Interestingly, it was found that spin squeezing is closely related to a key ingredient in quantum information theory 
and implies quantum entanglement [33l436l |. As there are various kinds of entanglement, a question arises: what kind 
of entanglement does spin squeezing imply? In Ref. [28], it has been found that, for a two-qubit symmetric state, 
spin squeezing is equivalent to its bipartite entanglement, i.e., spin squeezing implies bipartite entanglement and vice 
versa. Wang and Sanders [37|, presented a generalization of the results of Ref. (28| to the multiqubit case, where the 
authors showed that spin squeezing implies pairwise entanglement for arbitrary symmetric multiqubit states. More 
quantitatively, for spin states with a certain parity, if a state is spin squeezed according to the definition for K itagaw a 
and Ueda [15| , a quantitative relation is obtained between the spin squeezing parameter and the concurrence [42|-|45j , 
quantifying the entanglement of two half-spin particles [28l . . 

The close relation between the entanglement and spin squeezing enhances the importance of spin squeezing which 
motivate us to explore the role of Kerr medium [38j-|4l| and the nonlinear binomial state on the spin squeezing and 
entanglement. 

A binomial state is one of the important nonclassical states of light which has attracted much attention in the field 
of quantum optics in the last few decades, see for example |46h48| . This, in fact, is due to the importance of these 
states which have been experimentally produced and characterized. This state is regarded as one of the intermediate 
states of the coherent state [irj - lisl ]. It can be simply produced from the action of a single mode squeeze operator on 
the state | p, 0) [471 ]. where p is the Bargmann number. It also represents a linear combination from the number states 
with coefficients chosen such that the photon counting probability distribution is the binomial distribution with mean 
M\p\, where < \p\ < 1. 
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The scope of this communication is to employ the spin squeezing and quantum entropy to elucidate entanglement 
for two-atom system prepared initially in Bell-state interacting with single cavity field prepared initially in a binomial 
state. We introduce our Hamiltonian model and give exact expressions for the full wave function of the atomic Bell- 
state and field systems, shedding light on the important question of the relation between spin squeezing and quantum 
entropy behaviors. We also examine the evolution of the quasiprobability distribution Husimi Q-function for the 
model under consideration. With the help of the quasiprobability distribution, one gains further insight into the 
mechanism of the interaction of the model subsystems. In the language of quantum information theory a definition of 
spin squeezing is presented for this system. The utility of the definition is illustrated by examining spin squeezing from 
the point of view of quatum information theory for the present system. This analysis is applicable to any quantum 
tripartite-state subject to spin squeezing with appropriate initial conditions. 



II. MODEL SOLUTION AND THE REDUCED DENSITY OPERATOR 

Our system is as follows: A light field is initially in the binomial state (the state which interpolates between the 
nonlinear coherent and the number states ) 

oo 

\M0)) = \P,M) = J2 b n \n), (1) 

71 = 

where \n) is the eigenstate of number operator a) a = n, n\n) — n\n), while, the coefficients b n are given by 



1^77.(1 -H) M -", (2) 



where M > (i.e., M in general is any real positive number), < \p\ < 1, interacts with two identical two-level atoms 
are initially in the following Bell state 



|^ m (0))=7i|tt)+74|U>, (3) 



The binomial states have the properties 



M) 

\p,M) = {\0) |p|->0 (4) 
y) \p\ -> 0, M — > oo, M\p\ = a 2 . 

For this system, the atoms-field interaction is governed by the (N = 2) Tavis-Cummings model(TCM) 49]. The 
Tavis-Cummings model (TCM) ]49] describes the simplest fundamental interaction between a single mode of quantized 
electromagnetic field and a collection of N atoms under the usual two-level and rotating wave approximations when 
all of the atoms couple identically to the field. The two-atom (N = 2) TCM is governed by the Hamiltonian 

H = um + \u {d^ + + f( x , n) + A ]T + <ff a) , (5) 

where f{x, n) = xn(n — 1) + 2^/x ™ represents the nonlinear term with, 

fix, n) | n) = (xn(n - 1) + 2^/xn) | n) = fix, n)\n). (6) 



We denote by x the dispersive part of the third order susceptibility of the Kerr-like medium 38-41]. The parameter 
A represents the atoms-field coupling constant. The operators and ai l \ (i € {1, 2}) display a local SU(2) algebra 
for the i-th atom in the 2D supspace spanned by the ground (excited) state | J.),(| f)) and obey the commutation 

relations [a^ ,a^] = a[ l \ {i £ {1,2}), and a(a^) is bosonic annihilation (creation) operator for the single mode field 
of frequency cu. To make calculation to be more clear and simple, we put 

2 

k^2^P=Z; s = z,-,+, (7) 

3=1 
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the parameter k = 5, 1, 1 if s = z, — , +, respectively. 

The Hamiltonian (5), with the detuning parameter A = lu — to, takes the form 

H = H a + H int = u(n + J z )+AJ z + f(x, n) + X ( J_ a f + J+ a) , 



(8) 



Consider, at t = 0, the two atoms are in the Bell state, Eq. (3). The initial state of the system is a decoupled pure 
state, and the sate vector can be written as 



|^f(0)) = J2 6 ™ bi I tt. ») + 74 1 U, n)] , 



(9) 



n=0 



As the time goes, the evolution of the system in the interaction picture can be obtained by solving the Schrodinger 
equation 



ij t \ip AF (t)) =H\i> AF (t)), 

and using the condition Eq. (9) to obtain the time-dependent wave function in the form 

\^ AF {t)) = \U) \1) + \R) |2> + |S) |3) + |T) |4), 

where 



\U) = J2b n 



71=0 



A n (t)\n) + H n - 2 (t)\n - 2) 



(10) 



(11) 



(12) 



\R) = E 6 « 



n=0 



B„+i(t)|n+l) + G„_i(t)|n-l) 



= \S), 



and 



n=0 



D n+2 (t)\n + 2) +E n {t)\n) 



where we used the notations 

|l) = |tt), |2) = |tl> |3) = Ut) |4) = |U). 
the group of the complex amplitudes A n (t), D n+ 2(t), B n+ i(t) and C n +i(t) are given, respectively by 

1 2 

Mt) = ^jr ]T D k n+2 (0)[4 - 2Y\ + (a 2 + a 3 ) Vk + a 2 a^\ 
ll 11 2 '-^ 



(13) 



(14) 



(15) 



(16) 



fc=0 



1 2 

B n+ i(t) = -_£l>* +2 (0)(,fe + a 3 )e i ''* t = C n+1 (t), 



(17) 



with 



D n+2 {t)=Y. Dk ^y vkt i 



k=0 



(18) 



ai = A + /( X ,n), o 2 = /(x,n+l), a 3 = -A + /( X ,n + 2) 



(19) 



fi = Av/n+T, T 2 = AVn + 2, 



(20) 
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where 



with 



and 



3fe - ' xcos" 1 



X l 2 
IT + 3 



X 2 - 3X 2 cos(0 k ) 



9XiX 2 - 2Xl - 27X 3 



2(X 1 2 -3X 2 )5 



2kir\ 



k = 0,1,2, 



Xi = ai + a 2 + a 3 , X 2 = a 2 a 3 + a x (a 2 + a 3 ) - 2(T 2 + r|), X 3 = aia 2 Q!3 - 2 (a 1 r 2 . + a^r 2 ), 
where the complex coefficients D^ +2 (0), k = 0, 1, 2 are given by 

£>n+2(0) = 27lFlF2 ,fc,r,« = 0,1,2; fc ^ r ^ s 

Vkr Vks 

If in Eqs. (16-24), we let 71 — > 74 and 

ai = -A + /(x,n), a 2 = /(x,n-l), a 3 = A + /(x,n - 2), 



I\ = Ai/n, T 2 = AVn-T, 



(21) 

(22) 

(23) 

(24) 

(25) 
(26) 



and replacing the group of complex amplitudes A n (t), D n+2 (t), B n+ \{t) and C n+ \(t) by the other group of complex 
amplitudes E n (t), F n _i(t), G n -i(t), H n _ 2 (t), we obtain easily the last group, respectively. 
The reduced density operator of the subsystem is given by 



P{t)A{F) = T*F{A)P{t)AF = Tr F{A) \lp AF (t))(lpAF(t)\. 

Then the reduced atomic density operator in matrix form is given by 



PA 



where the elements p A r , s,r G {1,2,3,4} are given by 

pf = (U\U) = £{m 2 (k| 2 + \H n _ 2 \ 2 ) + b n b* n+2 A n H* n + b n b* n _ 2 A* n _ 2 H n _ 2 



(27) 



(pf 


Pa 2 


pf 


pf\ 


Pa 1 


P 2 2 


P 2 f 


Pf 


Pa 1 


Pa 2 


Pf 


Pf 


\pf 


pf 


Pf 


Pf) 



(28) 



(29) 



pf = (R\R) = f^\\b n \ 2 (\B n+1 \ 2 + |G„-i| 2 ) + b n b* n+2 B n+1 G* n+1 + &„6;_ 2 B;_ 1 G„_ 1 } = (S\S) = pf, (30) 

n=0 l- \ / > 

pf = (T\T) = J2\\b n \ 2 (\D n+2 \ 2 + \E n \A + b n b* n+2 D n+2 E* n+2 + b n b* n _ 2 D* n E n \, (31) 

n=0 l- \ / > 

00 , 

Pf = (R\U) = ^\b n b* n _ l A n B* n + b n b* n _ 3 B* n _ 2 H n _ 2 

n=0 ^ 

+b n b* n+1 A n G* n + 6„6:_ 1 G:_ 2J ff n _ 2 | = (p 2 iT = (S\U) = pf, (32) 
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n=0 



p\ A = (T\U) = V b n b* n _ 2 A n D* + b n b* n _ A D* n _ 2 H n - 2 



b n \ 2 A n E* n + b n b* n _ 2 E*_ 2 H n . 2 = ( P 4i)*, (33) 



pf = (S\R) =J2\ K\ 2 \B n+1 \ 2 + b n b* n _ 2 C* n _ x G n . x 

n=0 



+b n b* n+2 B n+1 F; i+1 + \b n \ 2 \F n _ x \ 2 = (R\S) = p 32 , (34) 



ft = (T\R) = b n b* n _ x B n+x D* n+x + 6„6; 1 _3^; 1 _ 1 G„_ 1 

n=0 



+b n b* n+1 B n+1 E* n+1 + b n b* n _ x E*_ x G n . x = (p i2 )* = pi 4 , (35) 



III. SPIN SQUEEZING 

Spin squeezing phenomenon reflects the reduced quantum fluctuations in one of the field quadratures at the expense 
of the other corresponding stretched quadrature. In the literature, [H, [HI, EH IH, [Ho|, HU , there are several definitions 
of spin squeezing and which one is the best is still an unsolved issue. Squeezing or reduction of quantum fluctuations, 
for arbitrary operators A and B which obey the commutation relation [A. B] — C, is the product of the uncertainties 
in determining their expectation values as follows [5l| : 

AAAB>±\{C)\, (36) 
where (AA) 2 = (A 2 ) - (A) 2 and (AS) 2 = (B 2 ) - (B) 2 . 

In this work spin squeezing parameters are based on angular momentum commutation relations. From the commu- 
tation relation [J x , J y ] = iJ z the uncertainty relation between different componenets of the angular momentum given 
by 

AJ x AJ y > ^\(J Z )\, (37) 



where 



J X = \{J+ + J-), (38) 



Jy = 1(J+ - J_), (.'!!)) 



where the operators J + , J_ and J z are given by (6). Without violating Heisenberg's uncertainty relation, it is possible 
to redistribute the uncertainty unevenly between J x and J y , so that a measurement of either J x or J y becomes more 
precise than the standard quantum limit y/\(J z )\/2. States with this property are called spin squeezed states in 
analogy with the squeezed states of a harmonic oscillator. Consequently, the two squeezing parameters can be written 
as 

Ft = (AJ X ) 2 l\(J,)\ = \(l- \Uj+ + J-))" \(J,)\) , (40) 
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F 2 = (AJ y f \\{J Z )\ = Ul- - J_)J> - |<J*)f), (41) 

If the parameter F% (F2) satisfies the condition F\ < (F2 < 0), the fluctuation in the component J x (J y ) is said to 
be squeezed. 

Using the wave function (11), we can easily compute the following time-dependent expectation values of the operators 
J_ , J+ and J z , respectively in the forms 

(J_)=2(^ 2 +^ 4 )^(J + )*, (42) 
W=P 1 A 1 -Pf, (43) 



IV. QUANTUM ENTROPY 

Quantum entropy, as a natural generalization of Boltzmann classical entropy, was proposed by von Neumann [52j . 
It has been applied, in particular, as a measure of quantum entanglement, quantum decoherence, quantum optical 
correlations, purity of states, quantum noise or accessible information in quantum measurement (the capacity of the 
quantum channel). Entropy is related to the density matrix, which provides a complete statistical description of the 
system. Since we have assumed that the two two-level atoms and the single-mode binomial field are initially in a 
disentangled pure state, the total entropy of the system is zero. In terms of Araki & Lieb inequality of the entropy (53j 

I S A (t) ~ S F (t) I < S AF (t) < I S A (t) + S F (t) I , (44) 

we can find that the reduced entropies of the two subsystems are identical, namely, SU(i) = S F (t). Here, the field 
entropy can be obtained by operating the atomic entropy. The quantum field entropy can be defined as follows |54| 

4 

S( PA ) = -J2 ni s) inn^, (45) 

s=l 

where IT^ , (s — 1,2,3,4) is the eigenvalue of the reduced density matrix, Eq. (27), and can be represented by the 
roots of forth order algebraic equation 

c + Cl n A + C2 u 2 A + c 3 n 3 A + n\ = o, (46) 

with the coefficients, c's, are given by 

C3 = -Pu - P22 ~ P33 ~ P44, (47) 

C2 = -|P4l| 2 - 2|/542| 2 - 2|pi 2 | 2 - |/023 | 2 + 2p 2 2(Pll + P44) + P22 + P11P44, (48) 
C\ = 2|pi 4 | 2 P22 + 2|p 2 4| 2 (Pll + P22) + 2|pi 2 | 2 (p22 + P44) - P22(Pll + P44) 
-2/911P22P44 ~ 3?(P23)(|P42| 2 + IP12H - < 2M{PilPl2P2i) + |p32| 2 (/Oll + Pte), (49) 
CO = |P4l| 2 |P23| 2 - P22P33|P4l| 2 - PllP33|P42| 2 - P44P22|P12| 2 - PllP44|P23| 2 - P33P44|Pl2| 2 
+P11P22P33P44 - P11P22|P24| 2 + (Pll|P24| 2 + P44|Pl2| 2 )5R(P23) + (^22 + P33 ~ 5R(P23) )^(P41 P12P24 ) , (50) 

where we used the notations 

P22=P33, Pl2=Pl3, P24=P34- (51) 

The eigenvalues 11^ , (s = 1, 2,3,4) are given, respectively by 

n y = El + (-l)'VH- 1 ; s=1)2 
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(a) _ U a + (-l) S+1 Vs+2 . 



ll A — 



s = 3, 4 



(53) 



where 



with 



C.3 



r.s = -y + (-) s f,Vs = ^/z 3 + (-iy Zi , 



(54) 



3c 2 8ci - 4c 2 c 3 + eg 
^3 = -2 c 2 + — - / , Zi = — 



(55) 



V. HUSIMI Q-FUNCTION 



For measuring the quantum state of the radiation field, balanced homodyning has become a well established method, 
it directly measures phase sensitive quadrature distributions. The homodyne measurement of an electromagnetic field 
gives all possible linear combinations of the field quadratures. The average of the random outcomes of the measurement 
is connected with the marginal distribution of any quasi-probability used in quantum optics. It has been shown from 
earlier studies |55T - l58j that the quasi-probability functions W-, (Husimi) Q-, and (Glauber-Sudershan) F-function, 
are important for the statistical description of a microscopic system and provide insight into the nonclassical features 
of the radiation fields. Therefore, we devote the present section to concentrate on one of these functions, that is 
the Husimi Q-function which has the nice property of being always positive and further advantage of being readily 
measurable by quantum tomographic techniques (59l. |60|. In fact, Husimi Q- function is not only a convenient tool 
to calculate the expectation values of anti-normally ordered products of operators, but also to give some insight into 
the mechanism of interaction for the model under consideration. The relation between the phase-space measurement; 
Husimi Q-function; and the classical information-theoretic entropy associated with quantum fields was introduced by 
Wehrl |6lj . However, on expanding the von Neumann quantum entropy in a power series of classical entropies, it 
was shown explicity [62[ that the first term of this expansion is the Wehrl entropy. Thus, Husimi Q-function can be 
related to the von Neumann quantum entropy in different approaches (6ll - [68| . 
The Husimi Q-function can be given in the form as (58l . l69l . PFoj 



Q{a) 



\<X\PF\ Ot) 



7T 



(56) 



where pf is the reduced density operator of the cavity field given by (27). The state | a) represents the well-known 
coherent state with amplitude a = X + iY. Inserting pp into Eq. (56), we can easily obtain the Husimi Q-function 
of the cavity field 



Q(a) = -Q (« I V) I 2 + I (« I T) | 2 ) 

7T 



(57) 



where 



n=0 



a a 
b n ^=A n {t) + b n+2 -=H n {t) 



(58) 



and 



r a* n „ . . a* n+2 
—=E n (t) + , 



D n+2 (t) 



(59) 



VI. DISCUSSION OF THE NUMERICAL RESULTS 



Using different sets of parameters for the initial state of the system we calculate numerically the quantum entropy 
Sa, spin squeezing parameters F% and F 2 and atomic population (a z ) as a reference function. All results are plotted 
as functions of the Rabi angle Xt. For each set of parameters four pictures are displayed. The pictures (a and b) show, 
respectively, squeezing parameters Fx and F 2 , while the pictures (c and d) show the quantum entropy Sa and the 
atomic population (a z ). For all our plots we set the Bell-state parameters 71 = -t= and 74 = 171. In figures I and 2 
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we have plotted the above mentioned functions with \p\ = 0.9, x/A = and various values of the detuning parameter 
A/ A. From these figures we can easily notice that, just after the onset of interaction these functions fluctuate for short 
period of time. This short period of revival is followed by a long period of collapse. The period of revival starts longer 
for one period of time with high amplitude to become wider with smaller amplitude as time goes. This is because 
the width and heights of the revivals belonging to the different series of eigenvalues are different. Furthermore, as 
we increase the detuning parameter A/ A from its value A/A = (resonance case), the overlap of revivals noticeably 
decreases with the increase in A/A. Also, the periods of oscillations within the revival decrease with the increase in 
the detuning parameter A/A. For the population, (er z ), the period of revival depends on the average number M\ p\ o f 
photons whereas the time of collapse depends on the dispersion, M|p|(l — \p\), in the photon number distribution [71j . 
Moreover, from these figures we can see that spin squeezing parameters Fi and F 2 experience collapse and revival 







FIG. 1: Spin squeezing parameters Fi (a), F2 (b), atomic entropy Sa (c), and atomic inversion (d) with |p| = 0.9, M = 50, 
X/A = and A/A = 



where atomic population exhibits collape and revival as time going. When we turn our attention to the role that spin 
squeezing parametres play to discover entanglement properties, we can realize that the behaviors of both squeezing 
parameter Fx and atomic entropy Sa are equivalent, i.e., entanglemet implies spin squeezing and vice versa. This 
can be understood as follows: quantum entropy Sa oscillates when F\ exhibits oscillations with same periods of time. 

Furthermore, the function Sa goes to its maximum when F\ shows oscillations around very small value (between 
0.45 and 0.5) of its maximum and when F2 shows collapse equal to its maximum, while Sa reaches its minimum 
when squeezing occurs. This behavior occurs periodically for both Sa and F±. This means that, on on-resonance 
atomic-system-field interaction, we can, with full success, understand entanglement dynamics from the dynamic of 
spin squeezing parameters F\ and F2 and vice versa. 

Let us now come to the case of off-resonance interaction between the atomic system and the cavity field. In this 
case the same general behavior (with periods shift to right when A/A = 10) is noticed. Additionally, the oscillations 
become more dense with reduced maximum of Sa corresponding to the increase of the oscillation interval of Fi 
(between 0.4 and 0.5 and become longer as A/A increase) and when some intervals of collapse begin to appear. 
The surprising and very interesting is the effect of the nonlinear medium individually and in the presence of the 
detuning parameter A/ A. To examine the effect of these mentioned parameters, we recall figures (3-5). These 
figures have been pictured by the setting of different values of the parameter \/X individually and in the presence 
of the detuning parameter A/ A. It is easy to see the change in figures shape where the standard behavior was 
changed completely. For a weak Kerr medium such that x/A = 0.5, our reference function, (er z ), shows behavior 
similar to the modified Jayned-Cumming model with Kerr medium (72l . [73| accompanied with reduced amplitude 
of oscillations. Furthermore, the population (cr z ) and spin squeezing parameter F% oscillate periodically with fixed 
periods are equivalent but F\ does not. A quick look at the squeezing parameter F\ one can realize easily that it 
oscillates rapidly. The oscillations overlap for a period of time (except for some instances) to become dense to show 



9 




50 100 






FIG. 2: The same as Fig. (1) but for A/A = 10 
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FIG. 3: The same as Fig. (1) but for x/A = 0.5 



periodically wave packets of Gaussian envelope with amplitude decreases as time develops. The more the Gaussian- 
packet-envelope amplitude decreases the more the entropy increases, i.e., stronger entanglement can be showed, see 
figure 3. This behavior becomes more clear when we consider the detuning parameter in our numerical computations, 
see figure 5. 

Furthermore, when the nonlinear medium becomes stronger, (a z ), F\ and show chaotic behavior with no indications 
of revivals or any other regular structure. This is accompanied with change in the entropy maximum from slow to 
rapid increase as x/A increases with time develops (see Fig. 4). This behavior is dominant without or with high 
values of A/ A. 
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FIG. 5: The same as Fig. (3) but for A/A = 5.0 



With the increase of \p\ and M, such that \p\ = 0.98, and M = 100, i.e.; on increasing the average number M\p\ 
of photon, the oscillations of squeezing parametres Fi, F2, (c z ) and quantum entropy Sa becomes more dense. 
This means that every two neighbor overlapping revivals start to overlap again. This is because with larger mean 
photon M these functions have bigger values with time evolution, which causes dense oscillations of the cavity 
field parametres. In other words, there are more revival series since more eigenvalues can be found in this model. 
However, the same behavior, we saw before for \p\ = 0.9, is seen again except for the envelope width becomes wider 
which resulted in fewer packets of oscillation appear in the same period of time, see figure 6 and 7. It is worth to note 
that each revival series of oscillations corresponds to a beat frequency. Moreover, on weak Kerr medium the relation 
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FIG. 6: The same as Fig. (1) but for \p\ = 0.98, M = 100 





FIG. 7: The same as Fig. (6) but for xA = 0.5 



between the atomic entropy and spin squeezing parametres seems more complicated. In this case all of them oscillate 
with no indications of revivals or any other regular structure where the Gaussian envelope completely disappears. 
The effect of strong Kerr medium individually and with the coexistence of small detuning exhibits behavior similar 
to that when \p\ = 0.9. 

At the end we are going to focus our attention on the representation of the field in phase space which provides 
some aspects of the field dynamics. Figure 8 shows mesh plots (left) and the corresponding contour plots (right) of 
the Husimi Q-function in the complex a-plane X — Re(a), Y = Im(a) for the Rabi angle Xt — ir/4 and different 
values of Kerr parameter x/A, while all other parameters are kept without change as in Fig. 1. From figure 8a, it 
is clear that the state of the field is a squeezed state, since the Husimi Q-function has different widths in the X 
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and Y directions. On the other hand, the squeezing is generated by the nonlinearity inherent to the system by the 
binomial state. This can be explained in terms of a superposition of different numbers of states of different phases, 
creating a deviation from the classical phase. It is well known that the squeezed states is a general class of the 
minimum- uncertainty states [74| . Bearing in mind that the nonlinearity of the binomial state yields squeezing of the 
quantum field and, as a result, entanglement between the two sub-systems is also produced [751 ] . It is of particular 
interest to see how the Husimi Q-function behaves once the Kerr medium is added. When x/A increases by a 
fraction value, we notice clearly the single blob become almost perfectly circular with radius \a\ « 7.5 rotates in the 
counterclockwise direction, see Figs. 8a, b, d and 8f. This behavior of the quantum field distribution is similar to 
that of the thermal state [74j which means that once the Kerr medium is added, makes it clear that rethermalization 
of the binomial field is indeed taking place. Now, it is perfectly sensible to ask, what is the situation if the Kerr 
Medium parameter is increased by an integer value? In this case, when the state evolves further, a multi-component 
structure develops as shown in Figs 8c, e and 8g, respectively. In these figures, the Husimi Q-function demonstrates 
that the quantum states obtained corresponding to Schrodinger cat states. Moreover, we see that the cat states have 
different number of components at different values of x/A. Different mechanisms demonstrated that, for the case of a 
radiation field propagating in a nonlinear medium, Schrodinger cat states are generated [66L l70l ] with different number 
of components at different times in the evolution |76| . It was shown that the splitting of the Husimi function, which is 
the signature of the formation of Schrodinger cat states, is related strongly to quantum entanglement (6rl Itm l77l - [79T | . 
To discuss the evolution of the Husimi Q-function in the case of resonance and fixed value of the Kerr parameter, 
i.e., x/A = 5.0 (strong Kerr medium), we have setted various values for the Rabi angle At in our computations, i.e., 
At = 0.0, 7r/6, 7t/4, 7t/3, 7r/2 and n. The results are displayed in figure 9. A collision of six blobs occurs gradually 
when At = 7r/6, which implies the rethermalisation of the quantum field. At the time evolution of At = tt, 7r/2,7r/3 
and 7r/4 the distribution of the Husimi Q-function splits into two, three and four blobs corresponding to Schrodinger 
cat states corresponding to different number of components at different times in the evolution. The center of blobs 
lies on a circle with radius \a\ centered at X = Y = 0.0. 
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FIG. 8: Husimi Q-function with \p\ = 0.9, M = 50, t = tt/4, A/A = and (a)x/A = 0.0, (b)x/A = 0.5, (c) x/A = 1.0, (d) 
X/A = 1.5, (e) x/A = 2.0, (f) X /A = 2.5 and (g) X /A = 5.0 
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FIG. 8: continued 



VII. CONCLUSION 

In Conclusion, we have shown that spin squeezing implies entanglement for quantum tripartite-state, where the 
subsystem includes the bipartite-state is identical. We have proved that spin squeezing parameters can be a convenient 
tool to give some insight into the mechanism of entanglement for the model under consideration. Moreover, a 
subsystem cavity field contains a nonlinear medium enhances noticeably the dynamics of entanglement specially 
when interacts with atomic subsystem off-resonantly. More clear insight into the relation between entanglement and 
the phase space distribution i.e., Husimi Q-function, is also illustrated. In this situation, the strong Kerr medium 
stimulates the creation of Schrodinger cat states which is necessary for the generation of entanglement. 
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FIG. 9: Husimi Q-function with \p\ = 0.9, M = 50, \ = 5.0, A/A = and (a) t = 0.0, (b)t = tt/6, (c) t = ir/4, (d) t = 7r/3, (e) 
i = tt/2 and (f) t = 7r 



[1] H. P. Yuen and J. H. Shapiro, IEEE Trans. Inf. Theory IT-24, 657 (1978); IT-25, 179 (1979); IT-26, 78 (1980); 

[2] C. M. Caves, Phys. Rev. D 23, 1693 (1981); J. Gea-Banacloche, and G. Leuchs, J. Mod. Opt. 34, 709 (1987). 

[3] X. G. Wang and B. C. Sanders, Phys. Rev. A 68, 033821 (2003). 

[4] U. V. Poulsen and K. Molmer, Phys. Rev. Lett. 87, 123601 (2001). 



-7.5 -5 -2.5 2.5 5 7.5 10 



FIG. 9: continued 

[5] I. Bouchoule and K. Molmer, Phys. Rev. A 65, 041803(R) (2002). 

[6] C. Genes, P. R. Berman and A. G. Rojo, Phys. Rev. A 68, 043809 (2003). 

[7] A. S. Sorensen and K. Molmer, Phys. Rev. Lett. 86, 4431 (2001). 

[8] L. Vernac, M. Pinard and E. Giacobino, Phys. Rev. A 62, 063812 (2000). 

[9] X. G. Wang, Opt. Commun. 200, 277 (2001). 

[10] X. G. Wang, Phys. Rev. A 331, 164 (2004). 

[11] D. Yan, X. G. Wang, L. J. Song and Z. G. Zong, Central Eur. J. Phys. 5(3), 367 (2007). 

[12] P. Zhou and J. S. Peng, Phys. Rev. A 72, 3331(1991). 

[13] A. Andre and M. D. Lukin, ibid. 65, 053819 (2002). 

[14] D. J. Wineland, J. J. Bollinger, W. M. Itano and D. J. Heinzen, Phys. Rev. A 50, 67 (1994). 

[15] M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993). 

[16] G. S. Agarwal and R. R. Puri, Phys. Rev. A 41, 3782(1990). 

[17] M.D. Lukin, S.F. Yelin, and M. Fleischhauer, Phys. Rev. Lett. 84, 4232 (2000). 

[18] A. Kuzmich, K. Molmer, and E.S. Polzik, Phys. Rev. Lett. 79, 4782 (1997) 

[19] J. Wesenberg and K. Molmer, Phys. Rev. A 65, 062304 (2002) 

[20] K. Helmerson and L. You, Phys. Rev. Lett. 87, 170402 (2001) 

[21] 6. E. Miistecaplioplu, M. Zhang, and L. You, Phys. Rev. A 66, 033611 (2002) 

[22] U. Poulsen and K. Molmer, Phys. Rev. A 64, 013616 (2001) 

[23] L. K. Thomsen, S. Mancini, and H.M. Wiseman, Phys. Rev. A 65, 061801 (2002) 

[24] T. Gasenzer, D.C. Roberts, and K. Burnett, Phys. Rev. A 65, 021605 (2002). 

[25] J. K. Stockton, J.M. Geremia, A.C. Doherty, and H. Mabuchi, Phys. Rev. A 67, 022112 (2003). 

[26] L. Zhou, H.S. Song and C. Li, J. Opt. B: Quantum Semiclass. Opt. 4, 425 (2002). 



17 



X. Wang, A. Sorensen, and K. Molmer, Phys. Rev. A 64, 053815 (2001) 
D. Ulam-Orgikh and M. Kitagawa, Phys. Rev. A 64, 052106 (2001). 

A. Dantan, M. Pinard, V. Josse, N. Nayak and P. R. Breman, ibib. 67, 045801 (2003). 

B. Julsgaard, A. Kozhekin and E. S. Polzik, Nature (London) 413, 400 (2001). 

V. Meyer, M. A. Rowe, D. Keilpinski, C. A. Sackett, W. M. Itano, C. M. Monroe and D. J. Wineland, Phys. Rev. Lett. 
86, 5870 (2001). 

D. J. Wineland, J. J. Bollinger, W. M. Itano, F. L. Moore and D. J. Heinzen Phys. Rev. A 46, R6797 (1992). 

A. Sorensen, M. L. Duan, J. I. Cirac and P. Zoller, Nature (London) 409, 1044 (2001). 

A. Sorensen, Phys. Rev. A 65, 043610 (2002). 

D. W. Berry and B.C. Sanders, New J. Phys. 4, 8 (2002). 

M. G. Raymer, A.C. Funk, B.C. Sanders, and H. de Guise, |quant-ph/0210137 



X. G. Wang and B. C. Sanders, Phys. Rev. A 68, 012101 (2003). 
M. S. Ateto, Int. J. Quant. Inf. 5, 535(2007). 
M. S. Ateto, Int. J. Theor. Phys. 48, 545 (2009). 

M. S. Ateto, Int. J. Theor. Phys. 49, 276 (2010); DOI 10.1007/sl0773-009-0201-0. 
M. S. Ateto, Applied Mathematics & Information Sciences, 3, 41 (2009). 
S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022(1997). 
W. K. Wootters, Phys. Rev. Lett. 80, 2245(1998). 

X. Q. Xi, S. R. Hao, W. X. Chen and R. H. Yue, Chin. Phys. Lett. 19, 1044(2002). 

G. F. Zhang and J. Q. Liang, Chin. Phys. Lett. 20, 452 (2003). 
V. Buzek, J. Mod. Optics, 37, 303 (1990) 

D. Stoler, B. E. H. Saleh, and M. C. Teich, Optica Acta, 32, 345 (1985) 
A. Vidiella-Barranco, J. A. Roversi, Phys. Rev. A 50, 5233 (1994). 
M. Tavis and F. W. Cummings, Phys. Rev. 170(2), 379 (1968). 

H. Pu, W. Zhang and P. Meystre, Phys. Rev. Lett. 89, 090401 (2002). 

D. F. Walls and P. Zoller, Phys. Rev. Lett. 47, 709 (1981). 
J. von Neumann, Gotingger Nachr., 273 (1927) 

H. Araki and E. H. Lieb, Commun. Math. Phys. 18, 160 (1970). 

V. Vedral, M. B. Plenio, M. A. Rippin and P. L. Knoght, Phys. Rev. Lett. 78, 2275 (1997). 

E. Wigner, Phys. Rev. 40, 749 (1932). 

Z. Wigner, Phys. Chem. B19, 203 (1932). 

K. E. Cahill and R. J. Glauber, Phys. Rev. 177, 1882 (1969). 

M. Hillary, R. F. O'Conell, M. O. Scully and E. wigner, Phys. Rep. 106, 121 (1984). 
E. L. Bolda, S. M. Tan and D. Walls, Phys. Rev. A 57, 4686 (1998). 
S. Mancini and P. Tombesi, Europhys. Lett. 40, 352 (1997) 
A. Wehrl, Rep. Math. Phys. 16, 353 (1979) 

V. Perinovd, J. Krepelka, J. Perina, J. Luks and P. Szlachetkta, Opt. Acta 33, 15 (1986) 
Hong-yi Fan and Qin Guo, |quant-ph/0611206fr l (2006) 



C. P.-Campos, J. R. G.-Alonso, O. Castanos and R. L.-Pena, cond-mat.guant-gas/0910.3256vl (2009) 



18 

[65] Li-yun Hu and Hong-yi Fan, quant-ph/0911.0125vl, (2009) 

[66] A. Miranowicz,J. Bajer, M . R. B. Wahiddin and N. Imoto, J. Phys A: Math. Gen. 34, 3887 (2001). 

[67] A. Miranowicz, H. Matsueda and M. R. B. Wahiddin, J. Phys. A: Math. Gen. 33, 5159 (2000) 

[68] G. P. Beretta, J. Math. Phys. 25, 1507 (1984) 

[69] K. Husimi, Proc. Phys. Math. Soc. Japan 22, 264 (1940) 

[70] H. Fu and A. I. Solomom, J. Mod. Opt. 49, 259 (2002). 

[71] A. Joshi and R. R. Puri, J. Mod. Opt., 34(11), 1421 (1987). 

[72] E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 (1963). 

[73] A. Joshi and R. R. Puri, Phys. Rev. A 45, 5056 (1992). 

[74] J. Rogel-Salazar, S. Choi, G. H. C. New and K. Burnett, cond-mat. guant-soft /03 02066 v 1 (2003) 

[75] K. Burnett, S. Choi, M. Davis, J. A. Dunninghman, S. A. Morgan and M. Rusch, C. R. Acad. Sci. IV 2, 399 (2001). 
[76] D. Yao, Phys. Rev. A 55, 701 (1997). 

[77] J. A. Vaccaro and A. Orlowski, Phys. Rev. A 51, 4172 (1995) 

[78] A. Orlowski, H. Paul and G. Kastelewicz, Phys. Rev. A 52, 1621 (1995) 

[79] I. Jex and A. Orlowski, J. Mod. Opt. 41, 2301 (1994) 



